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ABSTRACT 

Aims. The intensity profiles of the C 2 Swan bands in cool DQ white dwarfs cannot be adequately fitted with models that otherwise 
succesfully reproduce spectral features of the molecule CH in these stars. Initial modelling showed that a two-component atmosphere 
in the style of a spot might be able to solve the problem. 

Methods. We photometrically observed the two cool DQ white dwarfs GJ1 1 17 and EGGR78 to search for variability caused by stellar 
spots. 

Results. We have not found any such variability, but we estimate the effects of hypothetical spots on lightcurves. We also estimate 
detection probabilities for spots in different configurations. Alternative explanations of the problem are needed and briefly discussed. 
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1. Introduction 

Magnetism in white dwarfs (WDs) is usually studied through 
Zeeman splitting of atomic spectral lines. This method has 
proven to be very useful, and with the sensitivity added by spec- 
tropolarimetry it is very effective in detecting magnetic fields. 
Spectra of cool DQ WDs do not show any atomic lines, there- 
fore an alternative method is required. Fortunately, the carbon 
molecules present in most of these stars are sensitive to mag- 
netic fields and they can be studied using spectropolarimetry. 

During the past several years we have used circular spec- 
tropolarimetry to observe known DQ WDs to look for polariza- 
tion signals from C2 and CH molecules that can be found in their 
atmospheres. We have used two different telescopes, the Nordic 
Optical Telescope (NOT) on the Canary Islands and the ESO 
VLT (Cerro Paranal, Chile). So far, we have obs erved 12 objects 
and fo und one of them (GJ841B) to be magnetic (Vornanen et alj 
2010). The rest of the stars do not show any polarization signal 
at the noise level of our observations, i.e. 0.5 % in Stokes V/I 
for the NOT observations and 0.2 % for the VLT observations. 

In I Vornanen et al.l (201 Oh we s howed that the model pre- 
sented in Berdvugi na et al.l d2005l l2007t) works very well for 
the CH absorption bands in GJ841B and were also successful in 
modelling the CH and C 2 blend at 430 nm visible in both G99- 
37 and GJ841B. But while we modelled the intensity profiles 
of the purely C2 Swan bands, we were unable to achieve good 
fits. Problems with fitting Sw an bands in DQ W Ds have been re- 
ported before (for example, Dufour et al. 2005 have shown that 
the Av = +2 bands are consistently too strong in their model 
fits). Our model shows the same problem (as seen in Fig.[T|), but 
in addition to that the Av = band appear s too weak. Hints of a 
too weak Av = band can also be seen in Duf our et alj ((2005). 
These problems persisted for all non-magnetic DQ WDs in our 
survey, as long as we used a single temperature model. 

We started searching for a solution to this problem by com- 
bining two models in the manner of a photosphere with a spot 
to obtain reasonable fits to the spectra. This led to somewhat 
strange results. The best -fitting models usually had a small spot 
with T = 8000 - 10000 K and a photosphere of T = 2000 K. 
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Fig. 1. Model fit from lBerdvugina et al.l(l2005l) applied to GJ893, 
one of the non-magnetic objects in our survey. Transition desig- 
nations are given above the spectrum. The discrepancy between 
the depths of different Swan bands is obvious. The temperature 
of the model is 6000K. 



This would suggest a cool disk or an envelope of gas or dust 
around the WD with a hole in it through which the atmosphere 
would be seen. Although dust disks have been fo und around hot 
WDs and the central objects of planetary nebulae (Bilfkova et al. 
1201 ll) . this does not sound like a very probable scenario for our 
cool objects. Instead, a cool spot has been found b efore on a 
weakly (70 kG) magnetic WD (iBrinkworth et al.ll2005l) . 

To exclude spots or envelopes with holes as the reason for the 
observed properties in these WDs, we started a photometric ob- 
serving programme with a small remotely controlled telescope 
on the Canary Islands. Since WDs are rotating just like any other 
star, the movement of the spot on the visible hemisphere would 
cause some photometric variability. We report the results of our 
photometric observations here and discuss what we can deduce 
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Table 1. Observing log 
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Fig. 2. Lightcurves for EGGR78 and GJ1117 together with 
nightly averages (displaced for clarity). 



from them. At the end of the paper we discuss other possible 
reasons for troubles with modelling. 



2. Observations 

For our observations we used the 35cm Schmidt-Cassegrain 
Celestron telescope attached to the side of the 60cm KVA 
(Kungliga Vetenskapsakademien) telescope in Observatorio del 
Roque de Los Muchachos (ORM) on the island of La Palma 
in the Canary Islands. The telescope is remotely operated by 
researchers at Tuorla Observatory of the University of Turku 
(Finland). The small telescope has filters for BVRI photometry 
and is equipped with an Apogee Santa Barbara 47P CCD cam- 
era. The larger main telescope is used f or linear polarizatio n ob- 
servations with the Dipol-polarimeter (Piir ola et al.l l2005). The 
two telescopes are pointed to the same target and can be used 
simultaneously, if desired. 

We took data for GJ1 1 17 during seven nights spread over ten 
nights and for EGGR78 during 17 nights spread over one month. 
We took exposures in both B and R bands on all nights. The 
dates of the observations and number of exposures for each night 
are listed in Table [T] (available electronically only). Exposure 
times were 180 seconds forB-band and 100 seconds for/? -band. 
We performed differential photometry comparing the objects to 
three comparison stars in the field of view. Since there is no in- 
formation on the comparison stars apart from magnitudes in a 
few different bands (not UBVR), we also had to compare the 
comparison stars to each other to find out if any of them are 
variable. None of them were found to be variable within the er- 
rors. Figure [2] shows the lightcurves for EGGR78 and GJ1117 
in fi-band together with nightly averages compared to another 
star in the field of view. Errors of the individual measurements 
are from photon noise and the error on nightly averages are stan- 
dard deviations of all measurements for the night in question. 
The magnitudes are arbitrary. 
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Fig. 3. Lightcurves for EGGR78 and GJ1 1 17 from the last night 
of observations (June 2nd and March 31st, 201 1, respectively). 



The KVA telescopes are mostly used for optical monitoring 
of blazars in conn ection with the MA GIC Telescopes also lo- 
cated at ORM (See lAlbertetalJl2006l for example). The stabil- 
ity limit of the telescope is 15 mmag. As can be seen from Fig. 
HI there is no real variability in EGGR78 or GJ1 1 17 in B or in R 
band (not shown here). 

Figure[3] (available electronically only) shows as an example 
a lightcurve from the last night of observations for each star. No 
intra-night variability is evident in either object within errors. 
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Fig. 4. Periodograms for GJ 11 17 (top) and EGGR78 (bottom) 
show no signs of periodic variability. 



The periodograms (Figs. |4| are also devoid of any signs of 
variability apart from the roughly 1 day period from the observ- 
ing times. 

3. Results 

Obtained lightcurves appear to be flat within the errors an d nei- 
ther inspection by e ye, nor a Lomb-Scargle periodogram dLombl 
ll976tlScar gle 1982) can detect any kind of variability in our data 
set. 

Since we did not find any signatures of variability in our 
photometry, we considered carefully to which kind of varia- 
tions we are sensitive to and which we might be missing. We 
have two sources of uncertainty here. We might miss variability 
weaker than 15 mmag because of the limiting accuracy of our 
small telescope, and the timings of our observations also limit 
our sensitivity to certain rotational periods. 

To estimate the amplitude variations we could expect in the 
lightcurves, we used the model discussed earlier to calculate a 
lightcurve when a spot of higher C2 concentration and a tem- 
perature different from the atmospheric temperature crosses the 
visible hemisphere of the white dwarf. For simplicity, we as- 
sumed the spot to be a circular, flat disc. The angle between the 
normal of this disc and the line of sight is given by 



cos or = cos(/?)sin(/)cos(0) + sin(yS)cos(/), 



(1) 



where B is the spot latitude, i is the inclination of the rotation 
axis from the line of sig ht, and <j> is the rotational phase (<f>=0 to- 
wards the line of sight) (Bri nkworth et alj |2005). Conveniently, 
this is also the size of the spot projected onto the plane of the 
sky. We also ignored limb -darkening effects for simplicity. The 
model produces normalised intensity spectra, therefore we had 
to transform our observed lightcurve from magnitudes to fluxes 
and then normalise them to the average value of the flux to ob- 
tain comparable results from observations and modelling. 
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Fig. 5. Spot models: (I) homogenous atmosphere of 8000 K, (II) 
spot of 8000 K in an atmosphere of 2000 K, and (III) spot of 
6000 K in an atmosphere of 8000 K. 



We started by calculating the lightcurve from the 
spot/atmosphere combinations suggested by fitting the observed 
Stokes I spectra of the two WDs. The temperatures in these two 
models are 2000 K for the atmosphere and 8000 K and 10000 
K for the spots in EGGR78 and GJ1 1 17, respectively, with spot 
sizes of 1 % and 10 %. In both stars, modelling suggests that 
there is a low concentration of C2 in the atmosphere in general 
and the spot has a much higher amount of C2. Figure [5] shows 
three examples of model spectra: (I) homogenous atmosphere at 
8000 K, (II) spot of 8000 K in an atmosphere of 2000 K, and 
(III) a more realistic 6000 K spot in an atmosphere of 8000 K. 
The differences between the spectra are very small, especially in 
cases I and III. 

Since we do not know the rotational period of the white 
dwarfs, we used the time span of the photometric observations 
as the periods for illustration purposes. Results are shown in Fig. 

m 

It is evident from Fig. [6] that the spot model that fits the 
spectrum of EGGR78 will not produce a sufficiently deep dip 
in the lightcurve for us to detect it. It would take a size of nearly 
10 % of the visible hemisphere for the spot to be detectable. 
This is because the C2 Swan bands on EGGR78 are very weak 
and the filter sensitivity makes the variability of the absorption 
bands even harder to detect, since the strongest band between 
480 and 520 nm is not included in the B filter. Variability on 
GJ1 1 17 would clearly have been found if it were present. 

The model lightcurves in Fig.|6]have the most ideal spot con- 
figuration for detection: namely, a spot on the equator of the 
star. When the spot is moved towards the pole, the drop in the 
lightcurve becomes shallower as the projected area of the spot, 
and its variation, becomes smaller. The same happens if the ro- 
tational axis is tilted from the plane of the sky towards the ob- 
server. This also changes the duration of the dip because the spot 
spends a different amount of time behind the star. Finally, when 
the spot latitude B is higher than the inclination of the rotational 
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Fig. 6. Observed lightcurves (diamonds) of EGGR78 (top) and 
GJ1 1 17 (bottom) together with model lightcurves (solid line). 



axis i, the lightcurve assumes a true sine-form because the spot 
remains visible at all times but its projected size varies. 

To estimate our chances of detecting a spot depending on 
the system geometry, we calculated a lightcurve for all angles of 
spot latitude (-90° . . . 90°) and inclination (0° . . . 90°) in one de- 
gree intervals. If the difference between maximum and minimum 
values of brightness on the lightcurve was larger than 0.03, we 
counted it as detected. We also gave weights to the chances of 
detecting a spot for each latitude/inclination combination based 
on cos(yS - to account for the fact that the points appear to be 
more densely distributed away from the line of sight towards the 
edges of the stellar disk. In this way we estimated the chances of 
missing a spot on both stars. Using a spot size of 10 % we cal- 
culated that we would have had about a 40 % chance of missing 
a spot in both our WDs using either model II or III. 

There is, of course, the chance that we did not see variations 
in our photometric observations because the spot was behind the 
star during the observations. In this case, we might see some 
C2 that has spread out over the entire star, and missed the spot 
that is needed to model the spectra taken a year earlier. With 
the correct values of temperature and carbon abundance the spot 
can just alter the shape of the molecular absorption bands with- 
out changing their depths. The variations would be very difficult 
to pick up for the untrained eye by just looking at the spectra. 
Moreover, no one has tried monitoring these kinds of stars spec- 
troscopically for a longer period to see if they are variable or 
not. 

We should also consider the effect of a very long rotational 
period because our photometric observations as well as the spec- 
tropolarimetric observations might have been taken in the deep- 
est part of the brightness drop. The normalisation level of the 
model lightcurve is determined by our spectropolarimetric ob- 
servations, which were taken a year before the photometric ob- 
servations, therefore a long rotational period could mean that 
all our observations have been performed during a minimum. 
For example, a 3000 -day period would mean that both our spec- 



tropolarimetric and photometric observations could have been 
made in a low state. 



4. Conclusions 

We monitored two cool DQ WDs to search for photometric vari- 
ability as a signature of a stellar spot. We did not find any such 
variability. The accuracy and timing of our observations do not 
allow us to detect very long rotational periods (of about a year or 
longer), very small spots, or some spot configurations. Previous 
studies of WDs have come to the conclusion t hat the rotation 
perio ds usually last from hours to decades ( Charpine t et alj 
120091 and references therein). Some magnetic WDs have in- 
deed been found to be slow rotators w i th periods of about 
100 y ears dSchmidt & N orsworthvl Il99lt feerdvugi n & Piirolal 
1999). Since we did not detect any intra-night variability, the 
shorter end of the period distribution can be ruled out. Based on 
these arguments, we can assume that GJ1 1 17 and EGGR78 still 
might have very long rotational periods. 

Although we cannot rule out spots on these WDs com- 
pletely, we consider the values given by our model for the 
spot sizes and temperatures to be very unrealistic. But if the 
spot model is not the solution to the discrepancy in intensity in 
the carbon molecular bands, what better hypotheses could be 
found? We think the answer may lie in the oscillator strengths 
of carbon molecules that are embedded in our model. These 
parameters determine the strengths of the individual absorp- 
tion lines within the molecular bands. By modifying these values 
slightly, we hope to achieve consistent fits to all molecular bands 
without invoking a spot model. The current values were deter- 
mined in laboratory conditions, but still contain significant error 
limits. 

If modifying the oscillator strengths will not work, we have 
to go deeper into the physics of the problem and investigate in 
which way the structure, and therefore, the properties of a C2 
molecule is different in a white dwarf from conditions in the 
laboratory where the molecule was studied. Or maybe our treat- 
ment of the WD atmosphere is incorrect and we have to consider 
its properties more thoroughly. 
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